AD-A201  974 


fmrr' 


Deep-Sea  Research ,  Vol.  35,  No.  6,  pp.  973-983,  1988. 
Primed  in  Great  Britain. 


0198-0149/88  $3.00  +  0.00 
5)  1988  Pergamon  Press  pic 


h. 


i 

A  lee  wave  model  for  deep-sea  mudwave  activity 


Roger  D.  Flood* 

(Received  4  May  1987;  in  revised  form  19  November  1987;  accepted  25  November  1987) 


Abstract — Mudwaves,  with  wavelengths  up  to  6  km  and  heights  up  to  100  m,  are  commonly 
found  in  the  deep  sea  where  steady,  sediment-laden  currents  are  present;  their  internal  structure 
suggests  that  they  have  migrated  with  time.  Lee  waves  appear  to  be  generated  in  the  density 
gradient  above  the  sinusoidal  mudwave  topography;  the  near-bottom  flow  field  associated  with 
the  lee  waves  creates  a  cross-wave  asymmetry  in  bottom  current  velocity.  A  model  of  bottom  flow 
and  sedimentation  rate  for  a  transverse  mudwave  shows  that  preferential  deposition  occurs  on 
the  upstream  flanks  and  the  bed  forms  migrate  upstream.  The  flow  conditions  for  such  lee  waves 
are  common  in  the  deep  sea;  therefore  many  mudwaves  are  probably  active  under  present  flow 
conditions.  The  model  suggests  that  for  a  given  wave,  the  ratio  of  downstream-upstream 
sedimentation  rate  varies  primarily  with  flow  velocity.  Thus  changes  in  this  ratio,  determined  by. 
seismic  or  sampling  techniques,  might  be  used  to  determine  past  variations  in  flow  velocity.  (\  <  .. 


INTRODUCTION  AND  BACKGROUND 

Deep-sea  mudwaves  are  large-scale,  quasi-sinusoidal  bed  forms  found  in  many  ocean 
basins  where  deep  currents  play  an  important  role  in  sediment  transport  and  deposition 
(Fig.  1).  These  waves,  also  known  as  migrating  sediment  waves,  form  in  fine-grained 
sediments,  and  they  often  develop  on  the  flanks  of  large  sediment  drifts  (Fox  et  al. ,  1968, 
Hollister  et  al.,  1974),  along  continental  margins  (Jacobi  et  al.,  1975;  Lonsdale,  1983), 
and  near  deep-sea  channels  (Damuth,  1979;  Normark  et  al.,  1980).  The  heights  of  these 
waves  range  from  <10  to  >100  m,  and  their  spacings  range  from  <1  to  about  10  km.  The 
waves  can  be  followed  within  sub-bottom  profiles,  suggesting  that  the  waves  have  been 
present  for  very  long  time  intervals,  at  least  10’s  to  100’s  of  thousands  of  years. 

The  relationship  between  mudwaves  and  bottom  current  flow  has  been  poorly 
understood.  Detailed  surveys  of  mudwaves  related  to  contour-current  flow  often  indicate 
that  the  waves  are  typically  oriented  with  their  long  axes  at  a  low  angle  to  the  regional 
contours  (Hollister  et  al.,  1974;  Jacobi  et  al.,  1975).  The  waves  often  are  found  to 
migrate  upslope,  and,  because  the  waves  are  not  parallel  to  the  bottom  current  direction 
(which  is  presumed  parallel  to  the  regional  contours),  upeurrent.  In  some  instances,  such 
as  in  the  central  Argentine  Basin  (Flood  and  Shor,  1988),  the  mudwave  crests  appear  to 
be  perpendicular  to  the  bottom  current,  and  the  waves  migrate  upeurrent. 

Studies  of  mudwaves  on  the  Bahama  Outer  Ridge  (Flood,  1978)  suggested  that 
internal  waves  (lee  waves)  are  generated  in  the  weakly  stratified  near-bottom  flow  as  it 
passed  over  the  mudwave  (Fig.  2).  Studies  of  the  lee-wave  phenomena  by  Miles  (1968) 
and  Huppert  (1968),  among  others  (see  discussion  in  Turner,  1973),  suggest  that 
stratified  flow  over  an  isolated  ridge  will  create  strong  lee  waves  when  the  quantity 
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111 —  5  km  North  Flank,  Zapiola  Drift 

Fig.  1.  Examples  of  migrating  abyssal  mudwaves  from  the  Argentine  Basin.  Sub-bottom 
structure  on  these  3.5  kHz  records  suggest  that  the  bed  form  has  migrated  by  preferential 
deposition  of  sediments  on  one  side  of  the  wave,  perhaps  accompanied  by  erosion  on  the  other 
side.  Studies  of  mudwave  character  and  bottom  current  flow  directions  suggest  that  the  waves 
here  migrate  against  the  current,  that  is  that  preferential  deposition  occurs  on  the  upstream  side 
of  the  wave.  The  vertical  exaggeration  of  these  profiles  is  about  20:1.  The  waves  in  the  lower 
profile  appear  to  be  perpendicular  to  the  flow  whereas  those  in  the  upper  profile  are  at  an  angle  to 
the  flow  (Flood  and  Shor,  1988). 


OEEP-TOW  TEMPERATURE  SECTION  OVER  LARGE  MUD  WAVES 
BLAKE-BAHAMA  OUTER  RIDGE-tAREA  1) 
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Fig.  2.  Near-bottom  towed  temperature  profile  collected  across  a  mudwave  on  the  Bahama 
Outer  Ridge  (western  North  Atlantic)  near  28°17'N,  74°23'W.  (a)  Temperature  contoured  along 
Deep-Tow  fish  track,  (b)  Temperature  measured  at  Deep-Tow  fish.  During  the  time  of  this 
profile,  the  bottom  current  was  5  cm  s”1  (3.9  cm  s'1  perpendicular  to  the  wave),  and  k  =  NhIU 

was  0.4  (from  Flood,  1978). 
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k  =  N  h/U  (where  N  is  the  stability,  h  is  the  obstacle  height,  and  U  is  the  free-stream 
velocity  taken  normal  to  the  ridge)  is  near  unity  (Fig.  3).  For  flow  conditions  in  a 
mudwave  field  on  the  Bahama  Outer  Ridge,  k  values  range  from  0.4  to  1.5,  suggesting 
that  lee  waves  can  be  generated  under  present-day  flows.  Near-surface  currents  on  the 
downstream  flank  of  the  ridge  (where  the  streamlines  are  close  together)  are  stronger 
than  those  on  the  upstream  flank  (where  the  streamlines  are  farther  apart).  Flood  (1978) 
suggested  that  the  flow  pattern  could  lead  to  lower  sedimentation  rates  on  the  down¬ 
stream  mudwave  flank  and  higher  sedimentation  rates  on  the  upstream  side.  Such  a 
depositional  pattern  would  result  in  upstream  bed  form  migration. 

Such  a  lee-wave  model  of  mudwave  dynamics  can  be  quantified  in  the  case  where  a 
sinusoidal  mudwave  is  oriented  perpendicular  to  the  bottom  current  flow  to  predict 
variations  in  near-bottom  velocity,  bed  shear  stress,  and  sedimentation  rate  across  a  wave 
profile.  Queney  (1947,  1948)  presented  a  model  of  stratified  flow  over  a  sinusoidal  bed 
applied  to  the  problem  of  air  flow  over  mountain  ranges.  While  more  elaborate  models  of 
stratified  flows  over  topography  have  been  developed  since  then  (e.g.  Clark  and 
Peltier,  1984),  Queney’s  solution  is  of  particular  use  here  since  it  deals  with  flow  over 
sinusoidal  topography.  The  model  applies  for  small  deformations  of  the  flow  induced  by 
bottom  topography,  a  condition  met  by  mudwaves  as  the  wave  heights  are  almost  always 
<2%  of  their  wavelengths.  Since  the  model  applies  for  a  flow  which  is  perpendicular  to 
the  wave  crest,  it  will  not  apply  directly  to  rnudwaves  that  occur  at  an  angle  to  the  bottom 
flow  (e.g.  the  Bahama  Outer  Ridge);  however,  the  model  should  apply  to  areas  where 
the  waves  appear  to  be  perpendicular  to  the  bottom  flow  (e.g.  the  central  Argentine 
Basin). 


K  V 


Fig.  3.  Calculated  lee  wave  streamlines  created  by  a  stratified  flow  over  a  semicircular  ridge 
(Huppert,  1968).  Note  that  the  streamlines  are  more  closely  spaced  on  the  downstream  side  of 

the  ridge. 
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SEDIMENTATION  FROM  A  STRATIFIED  FLOW  OVER  SINUSOIDAL  TOPOGRAPHY 

For  a  rotating,  stratified  atmosphere,  the  stability  is 

N=  (g/pdp/dz)1'2,  (1) 

where  g  is  the  gravitational  acceleration,  p  is  the  density,  z  is  height. 

Queney  (1948)  presented  streamlines  for  stratified  flow  over  sinusoidal  topography  as 
a  function  of  k  (the  wavenumber  of  the  sinusoidal  profile,  k  =  2nlL,  where  L  is  the 
mudwave  wavelength),  ks  and  kf,  where  ks  =  NIU  and  kf  -  flU;  f  is  the  Coriolis 
parameter  and  U  is  the  mean  flow  velocity  of  the  stratified  layer. 

Solutions  were  presented  for  the  following  cases  (Fig.  4): 

Case  1,  ks  <  k:  the  vertical  variation  in  streamline  displacement  is  exponential  with  no 
phase  shift  over  the  wave. 

Case  2,  kf<  k  <  ks:  the  vertical  variation  in  streamline  displacement  is  sinusoidal 
with  a  phase  shift  over  the  wave. 

Case  3,  k  <  kf.  the  vertical  variation  in  streamline  displacement  is  exponential  with  no 
phase  shift  over  the  wave. 

For  mudwaves  in  the  deep  sea,  case  2  (lee-wave  formation)  leads  to  migrating 
sediment  waves.  For  any  given  wavelength,  lee  waves  form  from  Umin  -  ftk  to 
Umax  =  N/k.  Case  1  applies  where  the  flow  is  stronger  than  t/max  and  case  3  applies 
where  the  flow  is  weaker  than  £/„„ „.  The  flow  patterns  in  cases  1  and  3  lead  to  preferential 
deposition  in  the  trough. 

Following  Queney  (1948),  the  equations  for  the  streamline  displacement  C,  at  a  mean 
height  z  over  a  sinusoidal  seabed  topography  b  (for  case  2)  are  given  as  follows. 

(,  =  U  (B( po/p)1/2  cos(fct  +  k'z)  +  z),  (2) 

where  b  =  B  cos(fcx), 


ks<k  kf  <  k  <  ks  k  <  kf 

CASE  I  CASE  2  CASE  3 

Fig.  4.  Schematic  representations  of  streamlines  for  flow  over  sinusoidal  topography  in  the 
three  cases.  A  horizontal  offset  of  flow  lines  occurs  only  in  case  2.  In  cases  1  and  3  the  greatest 
streamline  spacing,  and  thus  the  lowest  velocity  and  highest  sedimentation  rates,  are  in  the 

trough  (Queney,  1947). 
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Here  U  is  the  mean  velocity,  B  is  the  mudwave  amplitude  (one-half  the  wave  height),  p  is 
the  mean  density  at  height  z,  po  is  the  density  at  z  =  0,  and  x  is  the  distance  along  the 
wave  profile. 

The  velocity  field  along  the  mean  flow  direction  can  be  determined  by  differentiating 
the  streamline  equation. 

u  =  U(  1  +  Bk'( po/p)1/2  sin(fct  +  k'z)).  (3) 

The  velocity  at  the  mudwave  surface,  z  —  0,  is 

u  =  U(  1  +  B  k'  sin(fct)).  (4) 

Since  k'  is  always  positive,  the  near-surface  velocity  field  reaches  a  maximum  value  at 
1/4  L  and  a  minimum  value  at  3/4  L.  The  vertical  gradient  in  u  is  small  near  2  =  0,  and 
zero  at  1/4  L  and  3/4  L,  and  the  equation  for  u  at  z  =  0  is  valid  for  z  up  to  several  10’s  of 
meters  at  1/4  L  and  3/4  L. 

The  equation  for  u  near  2  =  0  can  be  used  to  estimate  bed  shear  stress  across  the  wave 
profile.  The  bed  shear  stress  (r)  can  be  estimated  from  u  by  using  a  drag  coefficient  ( CD ) 
and  a  quadratic  stress  law, 

x  =  p  CD  u2.  (5) 

Drag  coefficients  have  been  determined  for  turbulent  Ekman  planetary  boundary 
layers,  where  u  is  the  free-stream  (geostrophic)  velocity  (e.g.  Csanady,  1967;  Weath¬ 
erly  and  Wimbush,  1980),  and  for  a  logarithmic  turbulent  boundary  layer,  where  u  is  the 
velocity  at  100  cm  above  the  bed  (e.g.  Sternberg,  1968).  Weatherly  and  Wimbush 
(1980)  determined  a  CD  value  of  0.002  for  a  deep-sea  turbulent  Ekman  boundary  layer 
whose  speed,  bottom  roughness  (z0),  and  stratification  are  similar  to  those  observed  near 
mudwaves  on  the  Blake-Bahama  Outer  Ridge.  Csanady  (1967)  suggested  CD  values  of 
about  0.001  for  atmospheric  boundary  layers  for  U/fz0  =  5  x  105,  conditions  similar  to 
those  of  Weatherly  and  Wimbush  (1980).  Drag  coefficients  determined  for  logarithmic 
turbulent  boundary  layers  range  from  about  0.002  to  0.004  (Sternberg,  1968).  Here  we 
assume  a  constant  CD  of  0.002  for  the  mudwave  profile  and  take  u  as  the  geostrophic 
velocity  above  the  Ekman  boundary  layer.  Weatherly  and  Martin  (1978)  suggested  that 
the  thickness  of  the  Ekman  boundary  layer  ( H)  can  be  represented  as  H  =  1.3  uj 
[f  (1  4-  A/2//2)1'4].  With  m*  =  0.36  cm  s-1  (corresponding  to  x  =  0.06  Pa),  N  =  5  x  10^* 
s'1  (a  typical  value  in  the  deep  sea)  and  at  45°  latitude,  H  is  20  m.  The  thickness  of  the 
logarithmic  turbulent  boundary  layer  is  0.10  -  0.15  H,  or  2-3  m  (Businger  and  Arya, 

1974).  For  simplicity,  velocities  calculated  at  z  =  0  m  are  considered  to  represent  the 
geostrophic  velocity  above  an  Ekman  boundary  layer  ( z  =  20  m).  This  approach  is 
justified  because  velocities  calculated  at  z  =  20  m  are  not  significantly  different  from 
those  calculated  at  z  =  0  m.  For  the  value  of  CD  used  here,  erosion  occurs  when 
u  =  17.1  cm  s-1.  Erosion  occurs  at  higher  velocities  for  smaller  values  of  CD  (for 
example,  erosion  occurs  at  24.6  cm  s-1  for  CD  =  0.001). 

McCave  and  Swift  (1976)  suggested  that  the  rate  of  fine-grained  sediment  deposition 
(/?,  g  cm-2  s“‘)  can  be  determined  from  the  sediment  concentration  in  the  water  ( Cw , 
g  cm3),  the  particle  settling  velocity  (Ws,  cm  s_1),  the  bed  shear  stress  (x),  and  the 
maximum  shear  stress  at  which  deposition  will  occur  (xc), 

R  =  CwWs(l-  x/xc)  P.  (6) 

$8  11  1^  062 


978 


R.  D.  Flood 


P  is  the  probability  that  a  sediment  particle  will  “stick”  when  deposited  (taken  as  1  here). 
Experimentally  determined  values  of  xc  range  from  0.04  to  0.08  Pa  (McCave  and  Swift, 
1976).  Here  we  take  an  average  value  of  0.06  Pa.  R  is  multiplied  by  3.16  x  1010  (seconds 
in  1000  y)  and  divided  by  1.5  (a  typical  near-surface  sediment  density)  to  convert  from 
g  cm2  s-1  to  geological  units  of  cm  y-1  x  10-3.  A  typical  value  of  Cw  for  regions  of  active 
sedimentation  in  the  deep  sea  is  1  x  10-7  g  cm-3  (100  pg  l-1;  Eittreim  et  al.,  1976; 
Biscaye  and  Eittreim,  1977).  McCave  (1985)  suggests  that  much  of  the  material  in 
suspension  at  the  HEBBLE  site  on  the  continental  margin  off  Nova  Scotia  falls  at  a 
velocity  of  5  x  10"3  cm  s_l.  As  t  increases  past  xc  sediment  ceases  to  deposit. 

We  assume  here  that  xc,  Cd,  Cw  do  not  vary  systematically  across  the  mudwave 
profile.  Near-bottom  studies  of  mudwaves  on  the  Bahama  Outer  Ridge  suggest  that 
bottom  roughness  is  similar  from  one  wave  flank  to  another  (Hollister  et  al.,  1974; 
Flood,  1978),  and  there  is  presently  little  evidence  to  indicate  that  the  concentration  or 
character  of  surficial  sediments  or  fine-grained  suspended  sediments  varies  systematically 
across  a  mudwave. 

The  maximum  sedimentation  rate  occurs  on  the  upstream  wave  flank  (3/4  L)  where 
the  flow  speed  is  the  lowest,  while  the  minimum  sedimentation  rate  occurs  on  the 
downstream  flank  (1/4  L)  where  the  flow  speed  is  the  highest.  By  forming  the  sedimen¬ 
tation  rate  ratio  (SRR)  between  the  minimum  (downstream)  sedimentation  rate  and  the 
maximum  (upstream)  sedimentation  rate,  both  Cw  and  W,  are  eliminated  from  the 
equation. 


xc-pCDU2{\  +  Bk'f 
xc-pCDU2{l-B  k’)2' 


(?) 


As  noted  above,  the  near-bottom  velocity  is  calculated  at  z  =  0  to  simplify  the 
expression.  If  CD  remains  constant  from  one  wave  flank  to  another,  the  sedimentation 
rate  ratio  (SRR)  will  depend  only  on  the  mean  flow  velocity  U  (k’  is  also  a  function  of  if). 


USE  OF  THE  MODEL  FOR  MUDWAVES  IN  THE  ARGENTINE  BASIN 

The  model  was  used  to  calculate  flow  and  sedimentation  rates  for  a  mudwave  with  a 
height  of  36  m  (amplitude  of  18  m),  a  wavelength  of  6000  m,  and  at  a  latitude  of  45°. 
These  parameters  approximate  those  found  in  the  field  of  apparently  transverse  mud¬ 
waves  in  the  central  Argentine  Basin  (Flood  and  Shor,  1988).  The  case  2  solution 
should  be  valid  for  a  mean  flow  of  9.8-47.7  cm  s-1 .  At  flow  speeds  slower  than  9.5  cm  s'1, 
sediments  should  accumulate  preferentially  in  the  wave  troughs.  Streamlines  calculated 
for  a  flow  of  13  cm  s-1  show  that  lee  waves  are  developed  above  the  mudwave  profile 
(Fig.  5)  and  k  =  0.13.  suggesting  that  the  interaction  is  weak.  The  near-bed  flow  speed 
(at  the  outer  edge  of  the  bottom  boundary  layer)  varies  from  14.3  cm  s_1  on  the 
downstream  wave  flank  to  11.7  cm  s-1  on  the  upstream  flank,  and  bed  shear  stress  varies 
from  about  0.041  to  0.027  Pa  (Fig.  6).  The  lowest  calculated  sedimentation  rate  is 
3.3  cm  y-1  x  10~3  on  the  downstream  wave  flank,  while  the  highest  rate  is  5.7  cm 
y'1  x  10-3  on  the  upstream  wave  flank  (Fig.  6).  The  sedimentation  rates  on  the  crest  and 
in  the  trough  are  equal  at  4.6  cm  y1  x  10"3. 

The  cross-wave  near-bed  velocity  and  sediment  rate  profiles  calculated  at  several  mean 
velocities  (Fig.  7)  show  that  the  sedimentation  rate  decreases  across  the  wave  profile  with 
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Streamline*  over  Slnueoldal  Topography  (I3cm/tec) 

>>>>>>  FLOW  DIRECTION  >  > 
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Fig.  5.  Calculated  streamlines  for  flow  over  a  mudwave  with  parameters  typical  of  those  in  the 
central  Argentine  Basin  (amplitude  =  18  m,  wavelength  =  6000  m,  latitude  =  45°)  at  a  mean 
velocity  of  13  cm  s'1.  The  upstream  shift  of  the  streamlines  with  height  above  bottom  suggests  lee 
wave  formation,  and  the  near-bottom  streamlines  are  most  closely  spaced  on  the  downstream 
flank  of  the  wave  and  most  widely  spaced  on  the  upstream  flank.  This  figure  has  a  vertical 

exaggeration  of  8. 
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Fig.  6.  Calculated  near-bottom  flow  speed  (c),  bed  shear  stress  (b),  and  sedimentation  rate  (a) 
for  a  flow  of  13  cm  s'1  over  sinusoidal  bottom  topography  (d,  mudwave  parameters  as  in  Fig.  5). 
The  highest  sedimentation  rates  occur  on  the  upstream  mudwave  flank  and  the  lowest  rates  occur 

on  the  downstream  flank. 


q  _  sedimentation  rate 


fo  _  bed  shear  stress 


g  _  flow  velocity 


increasing  flow  speed,  but  the  highest  sedimentation  rate  remains  on  the  upstream  wave 
flank.  For  a  mean  flow  greater  than  16  cm  s_1,  the  currents  on  the  downstream  flank  are 
strong  enough  to  erode,  and  the  zone  of  erosion  increases  with  mean  flow  speed  until 
sediments  can  deposit  only  on  the  upstream  flank.  At  higher  velocities,  erosion  occurs 
across  the  wave  profile.  For  these  parameters  the  sedimentation  rate  ratio  (SRR)  will 
vary  from  about  0.65  at  10  cm  s-1  to  0  at  16  cm  s_1  when  erosion  begins  to  occur  on  the 
downstream  wave  flank  (Fig.  8). 

EFFECT  OF  VARYING  WAVE  PARAMETERS 

The  changes  in  depositional  pattern  produced  by  changing  wave  dimensions  and 
latitude  can  be  most  readily  observed  by  plotting  SRR  for  waves  of  different  dimensions 
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mean  speeds  (cm/sac):  1  =  11  2  =  13  3=!5  4  =  17 
SRR  0.62  0.50  0.37  0.00 


Fig.  7.  Cross-wave  current  speeds  (a)  and  sedimentation  rates  (b)  calculated  for  several  mean 
flow  velocities  over  sinusoidal  topography  (c.  mudwave  parameters  as  in  Fig.  5).  Sedimentation 
rates  and  SRRs  decrease  as  flow  speed  increases.  Erosion  occurs  on  the  downstream  flank  at  a 
mean  speed  above  about  16  cm  s'1. 
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Fig.  8.  Plots  of  calculated  sedimentation  rate  ratio  (SRR)  vs  mean  flow  speed  for  mudwaves 
with  varying  amplitudes,  wavelengths  and  latitudes  (CD  =  0.002).  If  CD  is  taken  constant  at 
0.001,  the  SRR  vs  velocity  curves  are  similar  to  those  shown  in  Fig.  8;  however.  SRR  values  of  0 
occur  when  the  current  speed  is  about  24  cm  s  1 . 


(Fig.  8).  The  minimum  velocity  necessary  for  wave  migration  (t/mjn)  increases  with 
latitude,  suggesting  that  preferential  deposition  may  occur  in  high-latitude  wave  troughs 
even  at  moderate  flows.  Umin  also  increases  with  wavelength,  suggesting  that  mudwaves 
with  longer  wavelengths  will  become  inactive  at  lower  flow  speeds.  If  Umax  is  greater  than 
about  17  cm  s"1  (for  CD  =  0.002),  then  the  downstream  wave  flanks  can  be  erosional. 
For  N  =  5  x  KT4  s_1,  Umax  is  >17  cm  s-1  when  L  is  >2200  m.  Thus  preferential  erosion 
on  the  downstream  flank  is  more  likely  to  occur  on  mudwaves  with  longer  wavelengths. 

At  constant  velocity,  the  SRR  generally  decreases  with  increasing  wavelength  (com¬ 
pare  model  calculations  e,  b  and  a)  and  with  increasing  wave  amplitude  (compare  model 
calculations  d,  b  and  c).  Thus  mudwaves  with  longer  wavelengths  and  greater  heights 
tend  to  have  asymmetric  depositional  patterns  at  low  flow  speeds  (lower  SRRs),  whereas 
mudwaves  with  shorter  wavelengths  and  smaller  amplitudes  tend  to  be  more  symmetric 
even  at  higher  flow  speeds  (higher  SRRs) .  Mudwaves  at  higher  latitudes  are  somewhat 
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more  asymmetric  than  those  formed  at  lower  latitudes  (compare  model  calculations  f 
and  b). 

For  any  given  wave  parameters,  and  if  CD  is  the  same  on  both  mudwave  flanks  and 
with  time,  the  SRR  is  only  a  function  of  flow  velocity.  Since  the  SRR  is  also  the  ratio  of 
downstream  and  upstream  sediment  thickness  deposited  during  any  time  interval,  past 
values  of  the  SRR  can  be  determined  for  existing  sediment  waves  from  well-located 
piston  cores  or  high-resolution  seismic  profiles.  Such  information  may  allow  us  to 
calculate  bottom  water  flow  velocities  for  times  in  the  past  when  the  mudwaves  were 
active. 

The  assumption  of  constant  CD  across  a  mudwave  is  important  to  the  success  of  this 
sedimentation  model.  For  example,  if  CD  were  allowed  to  vary  across  the  wave,  then  CD 
on  the  downstream  flank  would  need  to  be  only  70-80%  of  the  value  on  the  upstream 
flank  in  order  to  have  equal  shear  stresses  (and  deposition  rates)  on  the  two  wave  flanks. 
However,  since  currents  on  the  downstream  wave  flank  are  predicted  to  be  stronger  than 
those  on  the  upstream  wave  flank,  current-created  roughness  is  probably  more  pro¬ 
nounced  on  the  downstream  flank,  and  thus  CD  on  the  downstream  flank  may  be  greater 
then  CD  on  the  upstream  flank. 

DISCUSSION  AND  CONCLUSIONS 

This  first-order  analysis  of  sedimentation  from  a  stratified  flow  over  a  sinusoidal 
topography  suggests  that  preferential  deposition  and  upstream  migration  are  natural 
effects  of  such  a  flow.  Preferential  deposition  occurs  on  the  upstream  mudwave  flank 
because  the  flow  pattern  creates  a  weaker  near-bottom  flow  there  (Fig.  9).  Less 
deposition  and  possibly  erosion  occur  on  the  downstream  flank  where  near-bottom  flow  is 
accelerated.  The  weak  near-bottom  density  stratification  of  the  flow  is  responsible  for  the 
development  of  such  forms  at  the  relatively  low  speeds  of  this  environment.  The  present 
model  is  notable  because  internal  waves  develop  on  the  density  gradient  that  is 
commonly  developed  at  the  base  of  the  water  column,  not  on  a  discrete  density  step  such 
as  might  exist  between  two  water  masses. 
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Fig.  9.  Summary  model  of  mudwave  growth  process.  The  lee  wave  flow  pattern  leads  to  lower 
bed  shear  stress  (higher  sedimentation  rates)  on  the  upstream  flank  and  higher  shear  stress  (lower 
sedimentation  rates)  on  the  downstream  flank. 
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The  model  simplifies  the  possible  flow  dynamics  of  mudwaves  considerably.  One 
important  simplification  is  that  only  the  cross-wave  flow  component  has  been  considered. 
Preliminary  investigations  of  the  along-wave  flow  component  (following  Queney,  1947) 
suggest  that  such  flows  can  be  significant  for  low  mean-flow  velocity.  However,  the  near¬ 
bed  flow  speed  is  <10%  greater  at  mean  speeds  >10-12  cm  s_1.  This  suggests  that  while 
there  may  be  significant  variations  in  local  flow  direction  across  the  wave,  the  cross-wave 
speed  variations,  and  the  SRR  plots,  are  similar  to  those  presented  here. 

Many  mudwaves  are  developed  at  an  angle  to  the  flow,  and  the  present  model  may  not 
be  directly  applicable  to  those  waves.  This  angle  may  be  the  result  of  a  long-term 
interaction  between  the  bottom  topography  and  the  flow.  Recent  studies  of  mudwaves  in 
the  Argentine  Basin  (Flood  and  Shor,  1988)  suggest  that  mudwave  orientations  can 
vary  regionally,  but  that  waves  on  south-facing  slopes  formed  by  eastward-flowing 
currents  have  similar  orientations.  Thus  effects  such  as  bottom  slope  or  the  earth’s 
rotation  may  be  important  in  controlling  wave  orientation.  Future  models  for  mudwave 
dynamics  will  need  to  include  cross-wave  flow  components  and  allow  for  topography  that 
is  developed  at  an  angle  to  the  flow. 

We  assume  here  that  the  same  lee  wave  flow  pattern  exists  for  a  long  period  of  time, 
and  thus  that  a  single  sedimentation  regime  applies  for  thousands  of  years.  Variability  in 
bottom  current  flow  patterns  may  be  due  to  tidal  currents,  eddies,  benthic  storms,  benthic 
fronts,  and  long-term  variations  in  bottom  water  flow  patterns.  All  of  these  components 
can  alter  the  relationship  between  the  wave  and  the  instantaneous  current  direction. 
Where  current  variabi'  n  is  small  (mean  flow  with  small  tides  and  few  benthic  storms),  a 
single  sedimentation  pattern  may  be  representative  of  long  periods.  However,  in  regions 
where  tidal  components  are  significant,  and  especially  where  they  dominate  the  mean 
flow,  the  simple  model  presented  here  may  not  apply.  The  passage  of  occasional  “benthic 
storms”  or  other  rare  flow  phenomena  also  will  affect  mudwave  activity  since  short-lived 
high-velocity  events  can  quickly  erode  sediments  deposited  over  much  longer  periods  of 
time.  Long-term  changes  in  bottom  water  flow  patterns,  velocity,  and  suspended 
sediment  load  and  character  also  must  affect  wave  development.  We  assume  that  bottom 
flow  and  sediment  characteristics,  including  CD,  tc,  Cw  and  Ws,  are  constant  across  a 
wave  profile.  The  validity  of  such  assumptions  needs  to  be  determined. 

If  the  predictions  of  this  or  a  more  refined  mudwave  growth  model  are  substantiated, 
then  past  SRRs  and  absolute  accumulation  rates  may  reflect  past  changes  in  mudwave 
activity  growth  which  can  be  related  to  changes  in  flow  rate  or  sediment  input. 
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